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20.  Abstract  (Continued) 

roll  rate,  induced  rolling  motion  with  a  priori  yawing  motion,  and  finally  the  combined  problem  of 
coupled  yawing  and  rolling  motion.  Critical  roil  rates  are  identified  and  solutions  valid  in  the  vicinity 
of  these  singular  conditions  are  generated.  Approximate  solutions  used  to  derive  stability  criteria 
are  compared  with  numerical  solutions  of  exact,  nonlinear  equations  of  motion.  Results  of  the 
perturbation  solutions  show  good  agreement  with  the  exact  calculations  for  moderate  angles  of 
attack. 
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NOMENCLATURE 


A  Missile  reference  area,  cm2. 

Cc  Total  missile  roll  moment  coefficient. 

Cp  Static  roll  driving  moment  coefficient. 

CCp  Linear  roll  damping  moment  coefficient. 

Cp  Induced  roll  moment  coefficient. 

Ns4 


(C;a  +  iCjj)  Complex  transverse  aerodynamic  moment  expressed  in  the  non-rolling  frame. 
Cm(y)  Roll-orientation-dependent  restoring  moment  coefficient. 

Cm  Restoring  moment  coefficient  derivative. 

Cm  Unmodulated  restoring  moment  coefficient  derivative. 

«n 


(Cm  +Cm  ) 
v  mq 


C  C 

'  mySi’  my6* 


Pitch  damping  coefficient  derivative. 

Coefficients  of  the  roll  orientation-dependent  static  moment  expansion. 


C 


m 


pa 


cn(t) 


CN 


764 


csnW 


Csm(T) 


Csma(T) 

Csm7S 


Linear  magnus  moment  coefficient  derivative. 

Roll  orientation-dependent  normal  force  coefficient. 

Unmodulated  normal  force  coefficient  derivative. 

Coefficients  of  the  roll  orientation-dependent  normal  force  expansion. 
Roll  orientation  dependent  side-force  coefficient. 

Roll  orientation-dependent  side-moment  coefficient. 

Roll  orientation-dependent  side-moment  coefficient  derivative. 

Linear,  roll  orientation-dependent  side-moment  coefficient  derivative. 


(Cy  +  iOj)  Complex  transverse  force  coefficient  expressed  in  the  non-rolling  frame, 
d  Missile  reference  diameter,  cm. 


H  CN  -  k72(C*  +C*  ). 

o  4 

Ix,  Iy  Missile  polar  and  transverse  moments  of  inertia,  Kgm-m2. 
K( ,  K2  Nutation  and  precession  modal  vectors. 

K,  ,  K->  Initial  values  of  the  nutation  and  precession  modal  vectors. 

lO  4 O 

ka  Missile  axial  radius  of  gyration. 
kt  Missile  transverse  radius  of  gyration, 
kj ,  k2  Nutation  and  precession  modal  amplitudes, 
m  Missile  mass,  kg. 


p  Roll  rate,  rad/sec. 

P  pd/V. 

(p,  q,  r)  Missile  angular  velocity  expressed  in  non-rolling  frame. 


pss  Steady -state  roll  rate. 

Dynamic  pressure,  l/2pV2. 


R  «„.■ 


R  R/(W[  _  w2). 
S  yk;2c;m(y). 


s  Nondimensional  arclength. 


T  Cl  +ka'2C*  . 


T0  Fast  time  scale. 

T2,T4  Slow  time  scales, 

t  Time,  sec. 


vi 


(o,  v,  w)  Missile  velocity  vector  expressed  in  the  non-rolling  frame. 


V  Magnitude  of  the  missile  velocity,  m/s ec. 
(x,  y,  z)  Inertial  frame  coordinates. 

(X,  Y,  Z,)  Missile-fixed  coordinate  frame. 

(X,  Y,  Z,)  Missile  non-rolling  coordinate  frame. 


GREEK  SYMBOLS 

<*T  Total  angle  of  attack. 

7  Cross  flow  orientation  relative  to  cruciform. 

7  u/V. 

S  Magnitude  of  the  complex  angle  of  attack, 
e  Smallness  parameter. 

9X ,  6 2  Nutation  and  precession  mode  slowly  varying  phase. 

X  Detuning  parameter. 

2  Complex  transverse  angular  velocity  expressed  in  the  non-rolling  frame. 
£  Complex  angle  of  attack  expressed  in  the  missile-fixed  frame. 

£  Complex  angle  of  attack  expressed  in  the  non-rolling  frame. 

p  Air  density  kgm/m  3 . 

0  Missile  roll  orientation. 

0O  Missile  initial  roll  orientation. 

0(s)  Total  phase  of  the  planar  motion  amplitude. 

Total  phase  of  the  jth  modal  amplitude,  kj. 

0  0j  +  0j 

cjj  Frequency  of  j*1'  modal  amplitude  kj. 

Wj  dcjj/dT2. 
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MATHEMATICAL  NOTATION 


(  )'  Denotes  differentiation  with  respect  to  non-dimensional  arclength. 

(  ‘  )  Denotes  differentiation  with  respect  to  time. 

(-)  Denotes  matrix. 

(  )  Denotes  quantities  expressed  in  the  non-rolling  frame. 

(  )  Denotes  the  complex  conjugate. 

|  |  Denotes  the  magnitude  of  a  complex  or  vector  quantity. 

1m  Imaginary  part. 

Re  Real  part. 
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CHAPTER  1 


s 


INTRODUCTION 


The  angular  motion  of  a  finned  missile  has  long  been  the  subject  of  extensive  experimental  and  analytical 
investigations.  Much  of  this  work  has  been  a  part  of  the  normal  aerodynamic  design  phases  of  both  guided  and 
unguided  missiles.  However,  all  too  often,  the  more  serious  and  ambitious  of  these  investigations  are  undertaken  as 
a  direct  result  of  unexpected  flight  failures,  at  the  expense  of  time  and  resources  already  invested  in  the  design.  In 
hindsight,  these  flight  instabilities  are  often  found  to  arise  from  fundamental  inadequacies  in  the  missile  design  or 
failure  on  the  part  of  the  designer  to  recognize  an  important  mechanism  of  aerodynamic  coupling. 

One  of  the  most  severe  types  of  flight  instabilities  encountered  in  the  free  flight  of  a  finned  missile  is  described 
by  large -amplitude  angular  motion  when  the  missile  rolling  velocity  is  in  the  vicinity  of  the  missile  natural  frequency 
(fundamental  resonance).  This  undesirable  behavior  has  been  shown  to  arise  as  a  result  of  the  aerodynamic  coupling 
which  exists  between  the  vehicle  axial  and  transverse  aerodynamics.  This  so-called  “roll-yaw  coupling”  has  been 
responsible  for  the  anomalous  behavior  of  several  contemporary  configurations  and,  in  some  instances,  has  resulted 
in  unacceptable  dispersion  levels  so  as  to  severely  degrade  weapon  effectiveness. 

In  spite  of  the  vast  amount  of  literature  available  on  the  subject,  a  complete  and  thorough  understanding  of 
the  phenomenon  has  successfully  eluded  the  flight  dynamicist  for  many  years.  While  there  have  been  significant 
contributions  made  in  recent  years1"7  in  the  form  of  partial  solutions,  the  problem  to  date,  in  three  degrees  of 
freedom,  remains  unsolved. 

While  previous  works  have  treated  the  possibility  of  missile  resonant  behavior  in  the  vicinity  of  fundamental 
resonance,  they  have  all  had  one  basic  assumption  in  common:  the  existence  of  a  small  configurational  asymmetry, 
which  acts  as  a  harmonic  forcing  function  and  is  the  cause  of  the  subsequent  resonant  behavior . 

If  one  considers  the  resonance  phenomenon  in  a  somewhat  more  general  light,  it  seems  plausible  that  the 
inherent  lack  of  complete  axial  symmetry  introduced  by  the  presence  of  the  fins  has  the  potential  for  introducing 
resonant  behavior.  Missile  transient  response,  in  the  presence  of  roll-orientation-dependent  aerodynamics,  at  a 
number  of  critical  rolling  velocities,  is  precisely  the  subject  of  this  report.  By  an  analysis  of  the  yawing  motion, 
with  the  rolling  motion  known  a  priori  and  vice  versa,  an  understanding  of  coupling  mechanisms  introduced  by  the 
roll -dependent  aerodynamics  may  be  gained. 

The  general  case  of  the  angular  motion  of  a  symmetric  missile  having  nonlinear  roll  orientation-dependent 
aerodynamics,  cubic  in  angle  of  attack,  and  constant  roll  rate  is  treated  via  the  method  of  multiple  scales.  Using 
the  resulting  first-approximate  solution,  missile  transient  response  is  investigated  and  three  critical  roll  rates  are 
identified  at  zero  spin,  resonance,  and  one-half  resonant  spin.  Plane  autonomous  system  theory  is  used  to  exhibit 
the  nature  of  the  motion  and  its  dependence  on  initial  conditions. 

The  induced  rolling  motion  of  a  cruciform  missile  with  a-priori  yawing  motion  is  then  studied,  using  a  gener¬ 
alized  version  of  the  method  of  multiple-time  scaling  to  obtain  the  conditions  for  roll  entrainment.  Solutions  are 
shown  to  be  singular  at  zero  spin,  resonance,  and  one-half  resonant  spinjand  additional  solutions  valid  in  the  vicinity 
of  these  roll  rates  are  generated. 

Finally,  the  full  nonlinear  problem  of  combined  pitch-yaw-roll  motion,  with  roll  orientation-dependent  axial 
and  transverse  aerodynamics,  is  studied,  with  no  a  priori  assumptions  regarding  the  motion.  Cri.xal  roll  rates  are 
identified  and  useful  stability  criteria  are  derived. 


CHAPTER  II 


EQUATIONS  OF  MOTION  AND  AERODYNAMIC  MOMENT  EXPANSIONS 

A.  EQUATIONS  OF  TRANSVERSE  ANGULAR  MOTION 

In  order  to  introduce  the  subject  and  gain  insight  into  the  problem,  the  equations  of  motion  of  a  rigid,  cruci¬ 
form,  free-flight  missile  having  both  aerodynamic  and  inertial  symmetry  are  presented.  In  the  work  that  follows, 
use  will  be  made  of  three  frames  of  reference:  an  inertial  frame,  a  missile-fixed  system,  and  the  “aeroballistic”  or 
nonrolling  system. 

The  (x,  y,  z)  coordinate  system  of  Figure  1  is  an  inertial  frame,  with  freestream  velocity  V  directed  along  the 
x-direction.  The  missile-fixed  frame  (X,  Y,  Z)  is  an  orthogonal  system,  rigidly  attached  to  the  missile,  with  the 
X-axis  along  the  longitudinal  axis  of  symmetry,  the  Y-axis  containing  a  fin,  and  the  Z-axis  completing  the  triad. 
It  is  convenient  to  express  the  equations  of  motion  in  terms  of  dynamical  quantities  defined  in  the  “aeroballistic” 
system,  (X,  Y,  Z),  which  is  allowed  to  pitch  and  yaw  with  the  missile  but  does  not  roll.  Transformation  of  vector 
quantities  between  the  missile-fixed  and  nonrolling  systems  is  then  completely  specified  by  a  rotation  about  the 
missile  X-axis  through  the  roll  angle  $  defined  by 


(2.1) 


where  p  is  the  missile  roll  rate  and  t,  time. 

Many  excellent  developments  of  the  equations  of  motion  of  a  symmetric  missile  can  be  found  in  the  literature 
and  for  the  sake  of  brevity  will  not  be  repeated  here.  Instead,  the  nomenclature  and  results  of  Reference  8  will  be 
used  in  writing  the  equations  of  angular  motion  in  the  nonrolling  coordinate  system  in  terms  of  the  translational 
(u,  v,  w)  and  rotational  (p,  q,  F)  velocity  components  expressed  in  that  frame.  From  Reference  8  the  exact 
equations  of  missile  angular  motion,  with  constrained  center  of  gravity  may  be  written* 

I'  “  iT?  =  C~  +  iC~  (2.2) 


?'  -  iP£  =  k”2(C~  +  iC~) 

where  the  complex  angle  of  attack,  £,  and  the  complex  angular  velocity,  £,  are  given  by 

t  =  2  + 

'  V 


(2.3) 


(2.4) 


£  = 


(q  +  ir)d 
V 


’Assumptions  of  constant  velocity  and  the  omission  of  gravity  arc  implicit  in  the  use  of  Equations  (2.2)  and  (2.3). 


(2.5) 


2 


Figure  1 .  -  Coordinate  Systems 


V  =  u/V 


and 


k"2  = 

Kt 


md2 


(2.6) 

(2.7) 

(2.8) 


In  compliance  with  the  methods  of  Reference  8,  the  independent  variable  t  has  been  replaced  by  the  nondimensional 
arclength,  s,  defined  by 


s 


(2.P) 


with  (  )'  denoting  differentiation  with  respect  to  s.  The  complex  moment  coefficient  in  the  nonrolling  frame 
(C;j  +  iC^)has  been  multiplied  by  the  density  factor  (pAd/2m) 
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(c~  +C~)  =  (C~  +  iC~ ) . 

'  m  n  *  '  in  n  / 


(2.10) 


Now  for  small  angles  of  attack,  the  quantity,  7,  is  essentially  constant  at  a  value  of  unity  since 

72  =  1  -  62  (2.11) 

where  6  represents  the  magnitude  of  f.  Therefore,  if  attention  is  restricted  to  small-amplitude  motion,  the  effect  of 
the  geometrical  nonlinearities  7  and  y'ly  may  be  approximated  by8 

7  *  1  (2.12) 

and 

7 7 7  ~  0  .  (2.13) 

With  the  elimination  of  these  geometrical  nonlinearities,  the  only  task  remaining  is  a  specification  of  the  functional 
form  of  the  aerodynamics. 


B.  aerodynamic  force  and  moment  expansions 

In  the  case  of  a  cruciform  missile,  the  static  restoring  moment  is  not  only  a  function  of  the  angle  of  attack 
but  is  also  dependent  on  the  orientation  of  the  plane  of  the  angle  of  attack  relative  to  the  cruciform,  as  shown 
in  Figure  2.  Experimentally,  this  modulation  of  the  static  restoring  moment  with  roll  angle  is  observed  to  be  quite 
small,  for  most  low-aspect-ratio  missile  configurations.9'11  Typical  variations  of  static  restoring  moment  with  roll 
orientation  are  shown  plotted  in  Figures  3  and  4. 


AFT  LOOKING  FORWARD 

Figure  2.  RoU-Qrientation-Dependeru  Aerodynamics 
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In  addition  to  static  moment  variations,  experiment  has  also  shown  the  existence  of  a  roll-orientation-dependent 
side  moment  and  its  corresponding  side  force.  This  “induced”  side  moment  is  caused  by  the  lack  of  complete  axial 
symmetry  introduced  by  the  presence  of  the  fins.  Thus,  when  the  cross-flow  vector  impinges  on  an  asymmetric 
body  cross  section,  at  specific  roll  orientation  angles  a  side-moment  arising  from  this  asymmetric  pressure  distribu¬ 
tion  results.  For  typical  low-aspect-ratio  configurations,  wind  tunnel  tests  have  shown  these  side  forces  and  moments 
to  be  reasonably  approximated  by  harmonic  functions  of  the  roll  orientation  angle,  7. 11  Typical  variations  of  side- 
force  and  side-moment  with  roll  orientation  are  shown  in  Figures  3  and  4. 

Given  some  rather  basic  assumptions  concerning  the  origin  of  the  fluid  forces  acting  on  a  missile  possessing 
n-gonal  symmetry,  Maple  and  Synge12  have  shown  that  simple  rotational  and  reflectional  symmetry  considerations 
can  sometimes  place  unexpected  restrictions  on  the  functional  form  of  the  missile  aerodynamic  force  and  moment 
system. 

The  basic  hypotheses  of  the  theory  are  twofold.  The  first  is  the  aerodynamic  hypothesis;13  i.e.,  the  aerody¬ 
namic  forces  and  moments  acting  on  the  missile  are  completely  specified  by  the  instantaneous  translational  and 
rotational  velocity  of  the  missile  relative  to  the  fluid.  The  second  is  the  assumption  that  these  forces  and  moments 
may  be  represented  by  a  power  series  expansion  in  the  transverse  translational  and  angular  velocity  components  of 
the  missile.  In  recent  years,  the  validity  of  the  aerodynamic  hypothesis  has  come  under  scrutiny,  and  in  some 
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Figure  3. -Typical  Restoring  Moment  and  Side  Moment  Variations 
with  Roll  Angle  (References  10  and  11) 
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Figure  4. -Typical  Normal  Force  and  Side  Force  Variations  with 
Roll  Angle  (References  10  and  1 1) 
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instances,  modifications  to  the  theory  have  been  proposed.  Considerations  of  such  mechanisms  as  nose-generated 
vorticity  and  its  interaction  with  the  aft  body  have  suggested  the  possible  inclusion  of  terms  that  would  reflect  the 
recent  history  of  the  flow  as  well  as  its  present  state.  Nevertheless,  over  the  years,  the  theory  has  essentially  remained 
intact,  providing  the  mathematical  framework  for  the  basic  aerodynamic  force  and  moment  expansion. 

Only  the  basic  results  of  the  theory  are  presented  herein.  A  development  of  the  static,  roll-orientation- 
dependent  force  and  moment  expansions  may  be  found  in  the  appendix.  For  a  general  treatment  of  the  subject,  the 
reader  is  referred  to  the  original  work  of  Maple  and  Synge,12  as  well  as  a  somewhat  abridged  version  by  Zaroodny.'4 

The  essence  of  the  method  is  embodied  in  the  observance  of  the  aerodynamic  force  and  moment  expansion 
under  coordinate  transformation.  The  physical  reasoning  of  these  operations  may  be  summarized  as: 

a.  The  aerodynamic  force  and  moment  expansion,  in  the  nonrolling  system,  must  remain  invariant  under  roll 
transformations  through  an  angle  2rr/N,  where  N  is  the  number  of  planes  of  rotational  symmetry  (fins). 

b.  The  aerodynamic  force  ard  moment  expansions  are  “antisymmetric"  under  a  reflection  of  the  coordinate 
system  about  a  plane  containing  a  fin  (cruciform  missile). 
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If  these  operation*  are  performed  for  a  cruciform  missile  and  attention  is  restricted  to  the  static-roll- 
orientation-dependent  terms  of  the  first  harmonic  in  47,  the  following  force  and  moment  contributions  (to  fifth 
order  in  6)  result: 


and 


-V’ V‘VM’*c»,^V4,t- 


(2.14) 


(Cffl+Cs)  =  i 


-C  +  C  62e~i4>  +  Cm  .fiV4-*  +  . . . 

ma(,  m>62  m>6< 


k- 


(2. IS) 


Equations  (2.14)  and  (2.15)  offer  a  considerable  amount  of  information  concerning  the  form  of  the  roll-orientation- 
dependent  aerodynamics  of  a  symmetric  cruciform  missile. 

First,  it  can  be  seen  that  the  terms  linear  in  angle  of  attack  are  free  of  any  7-dependence.  From  this  fact.it 
follows  that  the  eflects  of  roll  orientation  must  be  limited  to  terms  cubic  in  the  angle  of  attack  and  higher.  Terms 
of  lower  order  in  5  are  incapable  of  discerning  the  presence  of  the  fins.  Secondly,  if  Equations  (2.14)  and  (2.1 5)  are 
expanded  into  real  and  imaginary  parts,  it  is  seen  that  for  expansions  third  order  in  6  and  lower,  the  “in-plane”  and 
“out-of-plane”  components  due  to  roll  orientation  must  be  equal.  For  example,  the  in-plane  component  of 
Equation  (2.15)  is  given  by 


1  in-plane 


=  -Cm%5  +  CmWb*  Cos47  +  C"S846S  Cos4t  + 


(2.16) 


while  me  out-of-plane  component  may  be  written 


"out-of-plane 


-c 


mW 


6 3  Sin  47  -  Cm  ^65  Sin  47  + 


(2.17) 


It  is  easily  seen  that  descriptions  of  the  roll  orientation-dependent  static  moment,  up  to  third  order  in  5,  must  have 
the  side  moment  and  modulation  component  of  the  restoring  moment  equal  in  magnitude. 

In  the  study  of  roll  orientation-dependent  forces  and  moments  that  follows,  terms  up  to  third  order  in  6  and 
harmonic  in  47  will  be  retained. 

With  the  static-roll-orientation-dependent  forces  and  moments  defined,  inclusion  of  aerodynamic  damping 
(Cmq  +  Cmi),  and  a  linear  magnus  moment,  Cmpa,  results  in  the  total  aerodynamic  force  and  moment  system  to 
be  considered. 


(Cy  +  iCj)  = 

(<£+*:)  - 


-C*  +  T*  fi2e_i4'>  t 

CN  X/  * 


~iC*  +  iCm  ^V'4’’’  +  C*  f 

ma  m7  62  mpa 


$  -  iC  ?  +  Cl  £ 

md  mq 


(2.18) 

(2.19) 


Incorporation  of  Equations  (2.18)  and  (2.19)  into  the  equations  of  motion  with  the  assumptions  of  constant 
velocity,  small  geometrical  angles  (7  =  1),  and  the  omission  of  gravity  terms8  yields  a  second-order  differential 
equation  in  the  complex  angle  of  attack 


+  (H  -  iP)?’  -  (M  +  iPT)i  =  (R  +  iPS  )62e'i4l|  +  3C^  jV*4^! )'h\  (2.20) 


1 


where 


and 


H*V 


-2 


M  =  k:2c* 

'  m. 


"2  r* 


T  =  Ck  +  ka  Cm 


pa 


R  =  -k< ; 'cw 


(2.21) 

(2.22) 

(2.23) 

(2.24) 


S 


r 


=  C 


(2.25) 


The  axial  and  transverse  radii  of  gyration  are  denoted  by  ka  and  kt  respectively,  and  (”)  denotes  the  complex 
conjugate.  Note  that  the  right-hand  side  of  Equation  (2.20)  is  composed  of  both  lift  and  moment  terms  which  are 
nonlinear  in  angle  of  attack.  These  nonlinearities  are  partially  a  result  of  the  complex  exponential  e~i4y,  as  well  as 
the  nonlinear  angle  of  attack  dependence,  as  expressed  in  the  Maple-Synge  expansion. 


C.  DYNAMICS  AND  AERODYNAMICS  OF  ROLLING  MOTION 

The  equation  of  rolling  motion  of  a  symmetric  missile  is  also  taken  from  Reference  8  and  may  be  written 
simply  as 


lxp  =  q„AdCs  (2.26) 

where  Cs  is  the  total  aerodynamic  roll  moment  coefficient,  which,  for  the  purpose  of  the  present  analysis,  is  taken 
as  the  sum  of  a  linear  damping  and  a  static  aerodynamic  roll  moment,  C^. 

pd 

Cs  =  CCp  y  +  CCg  (2.27) 

Now  if  the  Maple-Synge  symmetry  arguments  are  applied  in  the  case  of  the  roll  aerodynamics  of  a  cruciform 
missile,  it  may  be  shown  that  the  static  roll  moment  terms  must  take  the  form  (see  the  appendix) 

oo  oo  oc 

C8  =r  C0sS2t  +  Y  S4N  y  62liBNeSin4N7  (2.28) 

S  c=o  n=1  e=o 

While  this  particular  expansion  is  based  on  tetragonal  symmetry,  it  does  not  require  missile  reflectional  symmetry 
and  hence  can  be  used  to  account  for  fin  misalignment. 

Now  if  Equation  (2.28)  is  expanded  for  the  lowest -order  terms  (8  =  0)  and  terms  quintic  in  angle  of  attack  and 
higher  are  omitted,  it  may  be  incorporated  into  Equation  (2.26)  to  yield  the  differential  equation  of  rolling  motion 
to  be  considered 


8 


(2.29) 


Ixp  =  q^Ad  CtQ  +  Ctp  +  Cv  ^464  Sin  47  +  . . .  . 

where  the  coefficients  in  the  formal  Mapie-Synge  expansion  have  been  replaced  with  the  more  familiar  Cy- 
(roll  driving),  and  C„  (induced  roll).  It  should  be  noted  that  while  Equation  (2.29)  contains  the  lowest-order 
terms  which  describe") fie  induced  roll  moment,  the  52,  54  dependence  of  the  roll  driving  moment  coefficient  has 
been  omitted,  since  these  terms  display  no  roll  orientation  dependence.  Hence  their  retention  was  not  considered 
essential  in  the  present  study  of  induced  rolling  motion. 

Typical  induced  rolling  moment  coefficient  behavior  with  roll  angle  is  shown  plotted  in  Figure  5  together  with 
yaw  moment  coefficient  data  taken  on  a  cruciform  missile  configuration  at  subsonic  Mach  numbers.  Also  shown 
plotted  are  the  analytical  approximations  to  the  data  which  resulted  from  a  least-squares  fit  with  the  Mapie-Synge 
functional  form. 
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Figure  5. -Typical  Induced  Roiling  Moment  Coefficient  Behavior  (Reference  1 1 ) 
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CHAPTER  III 


ANGULAR  MOTION  WITH  PRESCRIBED  ROLLING  MOTION 
A.  NONRESONANT  SOLUTIONS 

The  full  nonlinear  problem  of  pitch,  yaw,  and  rolling  motion,  in  the  presence  of  roll  orientation-dependent 
aerodynamics  is  the  ultimate  objective  of  this  analysis.  However,  it  will  be  instructive  in  what  follows  to  partially 
decouple  the  motion  by  considering  the  transverse  angular  motion  with  the  rolling  motion  prescribed  a  priori.  While 
this  is  not  representative  of  the  free-flight  condition,  it  hopefully  will  exhibit  the  effects  of  the  roll-dependent  forces 
and  moments  on  missile  stability  and  provide  insight  that  will  be  useful  in  the  solution  of  the  full  coupled  problem. 

In  this  chapter  a  solution  is  sought  to  the  differential  equations  of  yawing  motion,  in  the  presence  of  roll 
orientation-dependent  aerodynamics,  where  the  roll  rate  is,  at  most,  slowly  varying  in  comparison  with  the  frequency 
of  the  yawing  motion. 

An  inspection  of  the  nonlinear  roll -dependent  terms  on  the  right-hand  side  of  Equation  (2.20)  reveals  that  the 
imaginary  part  of  the  coefficient  R  +  iPS7  is  quite  small  in  comparison  with  R  for  roll  rates  characteristic  of  a  finned 
missile,  and  hence  may  be  approximated  by  the  real  part  only.  In  like  manner,  it  can  be  shown*  that  the  lift  term 
containing  CN^2  is  several  orders  of  magnitude  smaller  than  the  moment  term  and  will  therefore  be  neglected. 

In  the  development  of  the  approximate  solutions  which  follow  it  will  be  instructive  and  more  physically 
meaningful  to  change  the  independent  variable  of  Equation  (2.20)  from  arclength  to  time  which  is  easily  accom¬ 
plished  through  the  use  of  Equation  (2.9)  and  the  assumption  of  constant  velocity.  With  the  above  approximations 
in  mind.  Equation  (2.20)  may  be  rearranged  and  written  in  the  time  domain  as 

?  -  iP0?  -  u>q?  =  +  +  Xi$2e“l47? 


where 

P0  =  plx/ly 

W2  =  MV2/d2 
e2p,  =  -HV/d 

(3.1) 

e2p2  =  P0TV/d 
Xi  =  RV2/d2 

62  =  U 

7Z,  xs 

•The  term  t  /  £  is  proportional  to  the  frequency  of  the  yawing  motion  which  for  normal  epicyclic  motion  is  usually  on  the  order  of 
10"  7  rad/cal. 


k 


tl 


and  e  is  a  small  dimensionless  quantity  of  order  6  which  has  been  introduced  into  the  linear  damping  and 
magnus  terms  such  that  the  forcing  terms  to  the  unperturbed  differential  equation  are  of  order  63. 


In  order  to  obtain  an  approximate  solution  to  Equation  (25)  with  roll-dependent  aerodynamics  and  constant 
spin  via  the  method  of  multiple  scales,  two  independent  time  scales  T0  and  T2  are  introduced  such  that 


t  =  T0  +  e2T2  (3.2) 

Here  the  fast  time  scale  T0  is  associated  with  the  yawing  motion,  while  the  slow  time  scale  T2  is  associated  with 
quantities  changing  slowly  along  the  trajectory  such  as  dynamic  pressure,  aerodynamic  coefficients,  etc.  Thus, 
derivatives  with  respect  to  time  t,  are  transformed  to  partial  derivatives  according  to 


_d_  _  9  2  9 

dt  3T^  £  9T2 


(3.3) 


A  solution  to  Equation  (3.1)  is  sought  in  the  form  of  an  asymptotic  expansion  in  the  smallness  parameter  e  of 
the  form: 


1(0  -  eii(To>T2)  +  e3§"3(To>  T2)  +  ... 


(3.4) 


which  remains  uniformly  valid.  This  is  equivalent  to  the  condition  that  £3/£j  remains  bounded  for  all  T0  and  T2. 
The  introduction  of  the  smallness  parameter  e ,  into  the  equations  of  motion  and  the  assumed  form  of  the 
asymptotic  expansion  for  |(t)  may,  at  first,  seem  somewhat  arbitrary.  However,  a  successful  application  of  the 
method  of  multiple  scales  transforms  the  nonlinear  differential  Equation  (3.1)  to  a  sequence  of  linear  partial  dif¬ 
ferential  equations  in  the  £j(t). 


Substituting  Equation  (3.4)  into  Equation  (3.1)  while  making  use  of  Equation  (3.3)  and  equating  equal  powers 
of  e  one  obtains 


til 

9Tn 


iPn 


a l 

3Tn 


“  w0^1  -  0 


(3.5) 


3**3 


3*3 


—  iP0  = 


9Tn 


~292?i 

9T03T2 


9*! 


3*1 


+  iP0  —  +  ~  +  ^2*1  +  *i($i)  e 


3  i4 <p 


9T, 


(3.6) 


The  general  solution  to  Equation  (3.5)  is  the  usual  epicyclic  solution 

*l(T„,T2)  =  K,(T2)eiw'T<>  +  K2(T2)eiW2T<>  (3.7) 

where  and  cj2  are  the  nutation  and  precession  frequencies  respectively,  given  by 


wl,2  “  2  (P0  1  4w0 j 


(3.8) 


Having  obtained  a  solution  to  Equation  (3.5),  Equation  (3.7)  may  be  substituted  into  Equation  (3.6)  to  obtain 


where 


3Tn 


9?3  2~ 

-  iP°  aF-  w°*3  =  w‘e 


itJiTo  + 


j2T0 


+  Xiei40o|K?ei(4p_3t"‘)To  +  2Kfic2ei(4p’2w''^>To 

+K2ic1ei<4p_2^-w.)To  +  ic2K2ei(4p"2“''^)To 

+2K1K^ei(4p_<J>'2^)To  +  k3ei(4p‘3w2)To} 


(3.9) 


dK, 

Wj  =  i(P0  -  2cjj)  +  i^ijKjWj  +  i/r2Kj  (3.10) 


and 

dK2 

w2  =  i(P0-2w2)  +  ipjK2w2  +  4*2K2  (3.11) 


Implicit  in  Equation  (3.9)  is  the  assumption  of  constant  spin  such  that 

<t>  -  0O  +  pT0  (3.12) 

Now  particular  solutions  to  Equation  (3.9)  contain  secularterms  of  the  form  T0eia,jT<),  which,  as  T0  becomes 
large,  cause  £3/^  to  become  unbounded.  Therefore,  to  obtain  a  uniformly  valid  expansion,  the  Kj  are  chosen  such 
that  the  secular  terms  vanish,  i.e.,  wj  =0.  The  complex  modal  vectors,  Kj  may  be  represented  by  their  respective 
amplitude  and  phase  and  substituted  into  Equations  (3.10)  and  (3.11)  while  separating  real  and  imaginary  parts  to 
obtain  the  modal  damping  and  phase  differential  equations: 


dkj 

(p,wi  +  /a2) 

(3.13) 

dT2 

P0  ~  2wj 

dk2 

(Mi«2  +M2) 

(3.14) 

dT2 

P0  "  2w2 

d0t 

dT2 

d02 

=  dT^  =  0 

(3.15) 

where 

Kj  =  kjei8i  (3.16) 
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(w/VHdeg) 


With  the  secular  terms  removed.  Equation  (3.9)  may  now  be  solved  for  f3  to  obtain 


?3(T0.T2)  =  -Xi- 


K3e"'iTo 

-  P0vj  +  W5 


K,V"-T° 


3KfK2eu'jT° 


3K,K]ei"4Xo 


v\  -  P0p2  +  u>l  i>]  -  P0u3  +  <4  vl~  +  wo. 


(3.17) 


where 

t»l  =  4p  ~  3wj 
t>->  =  4p  -  3co-> 

(3.18) 

t>3  =  4p  ~  2tO|  -  to2 
f4  =  4p  -  Wj  -  2w2 


Thus  the  complete  solution  ?(Tq,  T2)  up  to  order  e3  may  be  written 


f(T0,T2)  =  eKieiw'T»  +  -  eh  J  1 5 - - 

vl  Pq^I  +  w0 


K^eil’2T°  3K?K2e‘‘'3T°  3K1K.?ei‘'4T° 

+ - - - +  - - — - -  +  - — - f  (3.19) 

v2  -  Po1^  +  «0  v\  ~  P0^3  +  w0  vl~  Po"4  +  w0  J 

Equations  (3.13)  through  (3.15)  indicate  that  to  a  first  approximation,  for  nonresonant  motion,  the  behavior 
of  the  modal  amplitudes  is  unaffected  by  the  presence  of  the  roll-orientation-dependent  aerodynamics*  Figure  6 
shows  a  comparison  between  solutions  generated  with  Equation  (3.19)  and  a  direct  numerical  solution  of 
Equation  (3.1). 


Figure  6.  -  Comparison  of  Approximate  and  Exact  Solutions:  Nonresonant  Spin 


•Identical  results  may  also  be  generated,  using  the  method  of  averaging.  (Reference  IS) 


fc 
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Missile  aerodynamic  and  inertial  characteristics  used  in  the  calculations  may  be  found  in  Table  1 . 

It  should  be  noted  that  Equation  (3.19)  becomes  singular  when  ^  =  oj12.  A  close  inspection  of  possible 
singular  roll  rates  reveals  that  Equation  (3.19)  ceases  to  be  valid  at  zero  spin,  p  =  and  approximately  one-half 
resonant  spin  p  =  (3cjj  +  w2 )/4.  At  these  roll  rates,  i/j  approaches  either  wj  or  Cl>2,  which  introduces  additional 
secular  terms  into  Equation  (3.9)  which  were  previously  not  taken  into  account. 


B.  RESONANT  SPIN 

Solutions  which  are  valid  at  resonant  spin  may  be  generated  by  a  consideration  of  Equation  (3.9)  with  i>]  =  coj 
to  obtain: 


— J  -  iP0  Jr  ~  3?3  =[wl  +  XiKfel40°]elu,lT° 

3Tq  3T0 

+w2ei^To  +  X]ei40o[K3ei(4w'‘3^)To 

+3KfR2el(2u;i"u'2)To  +  3K1K^ei<3wi'2uJ2)T°]  (3.20) 

With  the  elimination  of  secular  terms,  the  damping  and  phase  differential  equations  for  the  variation  of  the  modal 
amplitudes  take  the  form 


i(P0  ~  2coj )  -pp  +  iMiKjto,  +  iM2Ki  +  Xie‘40°Kj  =0 


dT 


(3.21) 


dK2 

i(P0  _  2co2)  +  ip]K2co2  +  i/t2^2  =  0 


(3.22) 


With  secular  terms  eliminated,  Equation  (3.20)  may  now  be  solved  for  ^(Tq,  T2)  which  is  valid  only  at  resonance 
(p  =  cjj)  to  yield 

K23ein'T° 


i3(T0’T2)  =  X,e‘40o 


Lnl  "  P0nl  +  w0  n2~  P0n2  +  w0  n3  _  P0n3  +  Wq 


3KjK2ein2T°  3K1K^ein3T° 


- 


2  _ 


(3.23) 


where 


n,  =  4w,  -  3w-> 


n2  =  2w,  -  u>2 


n2  =  3wj  -  2w2 

The  complete  solution  |(T0,  T2),  valid  at  resonant  spin,  may  now  be  written 


(3.24) 


i  u;.T, 


«T0.T2)  =  eKje'"'10  +  eK2e  2  0  +  e3Xie 


0  +  e3y,el40° 


K3ein>To 


nl  P0»l  +  w0 


3K,K2e  2  0  3K.K2e 
+ 


n2  Pon2  +  w0  n3  Pon3  +  w0 


(3-25) 
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Table  1 .  -  Missile  Aerodynamics  and  Configurational  Characteristics 


MACH  NUMBER 

0.78 

FLIGHT  ALTITUDE 

6000  km 

REFERENCE  DIAMETER,  d 

27.3  cm 

MASS,  m 

272.1  kgm 

POLAR  MOMENT  OF  INERTIA 

2.7  kgm  -  m2 

TRANSVERSE  MOMENT  OF  INERTIA 

54.5  kgm  -  m 

AERODYNAMICS 

CN  =  4.30 
‘a 

^mq  + 

Cma  =  -4.20 

Cm  =  5.00 
Pa 

c^v  =  9-‘ 

C?  =  0.20 

0 

CnV2  =  ~30'8 

Cc  .  =  10.6 

Cn  =  -2.0 
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Several  important  comments  regarding  the  distinction  between  Equations  (3.19)  and  (3.25)  should  be  made. 
While  Equation  (3.19)  is  valid  for  all  spin  except  p  =  0,  to,  12,  tO) ,  Equation  (3.25)  is  valid  at  p  =  Wj  only.  For  non¬ 
resonant  spin,  the  effect  of  the  roll-dependent  aerodynamics  was  not  discernible.  However,  at  resonance  (p  =  ojj), 
the  nutation  amplitude  may  be  strongly  affected  by  Xt  and  This  characteristic  of  sensitivity  of  the  stability  of 
the  motion  to  initial  conditions  gives  some  insight  into  the  mechanism  of  resonance-related  instabilities  whereby 
the  missile  stability  is  dependent  upon  the  amplitude  as  well  as  the  initial  phase  of  the  motion. 

The  validity  of  Equation  (3.25)  was  checked  by  comparison  with  direct  solutions  of  Equation  (3.1)  at  resonant 
spin  and  the  results  of  these  calculations  are  shown  in  Figures  7  through  10  where  two  identical  configurations  are 
released  at  slightly  different  initial  conditions.  Figure  8  is  seen  to  depict  stable  motion  while  the  motion  of  Figure  10 
shows  the  configuration  to  be  slightly  unstable. 


Additional  insight  into  the  sensitivity  of  the  motion  at  resonant  spin  to  variations  in  initial  conditions  may  be 
gained  through  a  decomposition  of  Equation  (3.21)  and  displaying  the  slowly  varying  dynamics  of  the  two- 
dimensional  system  in  an  amplitude-phase  plane.  Separating  Equation  (3.21)  into  real  and  imaginary  parts,  the 
nutation  mode  behavior  may  be  represented  by 


dkj 

dT7 


kj  Sin  0O 


(3.26) 


dfl, 

ir2 


k)  Cos  0O 


(3.27) 


Figure  7.  -  Complex  Angle  of  Attack  History  at  Fundamental  Resonance  with 
Nonlinear,  Roll-Dependent  Aerodynamics 
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Figure  10.  -  Missile  Transient  Response  at  Fundamental  Resonance  with 
Nonlinear,  Roll-Dependent  Aerodynamics 


where 

h  =  ^<t>0~eO  (3-28) 

Solutions  to  Equations  (3.26)  and  (3.27)  are  shown  as  trajectories  in  the  k\  -0O  plane  of  Figure  1 1 .  Initial  condi¬ 
tions  leading  to  stable  motion  are  seen  to  be  those  originating  below  the  dashed  line  denoted  as  the  separatrix.  It  is 
now  clear  that  the  principal  features  determining  the  stability  of  the  missile  at  resonance  are  the  amplitude  of  the 
nutation  component  kj ,  and  the  relative  phase,  0O. 


k^(rad)2 


CHAPTER  IV 

INDUCED  ROLLING  MOTION  WITH  PRESCRIBED  YAWING  MOTION 


In  the  previous  chapter,  attention  was  focused  on  the  nonlinear  yawing  motion,  while  the  rolling  motion  was 
constrained  at  constant  roll  rate.  While  this  can  hardly  be  considered  representative  of  free-fl ight  conditions, 
considerable  insight  was  gained  into  the  effect  of  the  roll-dependent  side  moment  on  missile  stability.  In  this 
chapter,  this  philosophy  is  applied  in  a  study  of  the  induced  rolling  motion  where  a  solution  to  the  nonlinear  roll 
dynamics  is  sought  while  the  yawing  motion  is  constrained  to  be  epicyclic. 


A.  NONRESONANT  SOLUTIONS 


The  roll  differential  equation.  Equation  (2.29),  which  was  previously  derived  in  Chapter  II  to  include  the  lowest  - 
order  roll-dependent  terms  consistent  with  the  Maple-Synge  expansion,  may  be  rewritten  in  slightly  different  form 
with  the  introduction  of  the  smallness  parameter  e.  whose  magnitude  is  of  the  order  of  the  total  angle  of  attack,  as 


where 


dp 

dt 


-> 

=  e~ 


°0  +  a2P 


e2o0  = 


q„Ad 
_ —  r 
lv  C«o 


(4.1) 


2  _  9„Ad  _ 


1  "  lv 


„Ad2 


VL  ~«p 


Approximate  solutions  to  Equation  (4.1)  are  sought  using  the  generalized  version  of  the  method  of  multiple  scales.5 
If  Equation  (4.1)  is  considered  in  the  absence  of  the  induced  roll  term,  the  solution  exhibits  the  usual  slowly  varying, 
exponential  behavior  of  the  roll  rate,  which  is  associated  with  the  slow  time  scale  T2  such  that 


d  _  2  9 
dT  "  e  TT2 


(4.2) 


Inclusion  of  the  induced  roll  term  introduces  a  nonlinear  forcing  function  whose  excitation  will  be  dependent  upon 
the  amplitude  and  frequency  of  the  assumed  missile  angular  motion,  ij(t).  If  the  familiar  epicyclic  motion  of  a 
symmetric  missile  is  assumed,  it  is  evident  that  the  resulting  solution  p(t)  should  contain  components  whose  fre¬ 
quencies  are  linear  combinations  of  the  epicyclic  frequencies  Wj  and  Wj.  In  view  of  this,  it  is  assumed  that  £(t)  is 
given  by 


f(t)  =  £  Kje‘  *1  (4.3) 

1  =  1 
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where 


d*j 


dt 


(4.4) 


and  the  asymptotic  expansion  for  p(t;e)  is  of  the  form 

p(t;  e)  =  p0(T2)  +  e2p2(i/t,,^2.0.T2)  +  •••  (4.5) 

Thus  time  derivatives  of  p(t)  will  be  converted  to  partial  derivatives  with  the  independent  variables  i//j,  <t>,  and  T2. 

Substitution  of  Equation  (4.5)  into  Equation  (4.1)  while  making  use  of  Equations  (4.2)  and  (4.4),  and 
equating  all  terms  of  order  e2  results  in  the  linear  partial  differential  equation 

aP2  dp2  3p2  r  4  i4(0,-0) 

+  "2  W2  P0  =  W3(T2>  [  16 

+K^e,4(*a'*)  +  4K2K2ei(3**ni_40) 

+6K2K2e2i(0'+,^20) 

+4K1^ei(J''+3^'40)]  (4.6) 


where 


dPo 

w3(T2>  =  _dT2  +  a2P0  +  °0 


(4.7) 


and 


K,  •  k/i 


(4.8) 


Particular  solutions  to  Equation  (4.6)  become  unbounded  as  ipj  and  <p  become  large.  Therefore.it  is  required  that  the 
secular  term  w3  vanish,  yielding  the  differential  equation 


dPo  _  x 

TfT  -  ao  a2po 


(4.9) 


for  the  slowly  varying  spin.  With  the  exclusion  of  the  secular  term,  a  particular  solution  to  Equation  (4.6)  may 
be  found 
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«•>  ■  4(S^t  c“4(J'  -*> 4  4isrk> Cos4tf  >  -« 


3o]k2k2 


4a,  k?k2 

+  - 7 - Cos  (3  $  i  +  $-,  -  40)  + 

4po  3a)j  w2  *  2pq  Wj  co2 

4a,k,k3 

+  W-UI  -3^2  Cos^i+3^-'40) 


Cos  2(0j  +  02  -  20) 


(4.10) 


Solving  Equation (4.7)  for  Po(T2),  the  complete  expansion  for  p(t),  with  induced  rolling  moment  in  the  presence  of 
epicyclic  yawing  motion  becomes 


P(t;e)  =  Po(°)e°212  +  -Jj£  (e0^2-!)  +  e2o, 


L4(Po“  wi) 


Cos  4(0,  -0) 


4 


4k]  k2 


+  77 - r  Cos  4(02  -  0)  +  - - - -  Cos  (30,  +  0-,  -  40) 

4(Pq  -  w2>  4pQ-3o;1-aj2 


3kjk? 


2po  ~  w1  -  w2 


Cos  2(0i  +  02  -  20)  + 


4k,  k 


1k2 


4Pq 


lo1  -  3to2 


Cos  (0 1  +  302  -  40) 


(4.11) 


Inspection  of  Equation  (4.1 1)  again  reveals  the  existence  of  discrete,  singular  roll  rates  at  which  the  solution  is  no 
longer  valid.  If  the  denominators  of  the  e2-tcrms  of  Equation  (4.1 1)  are  allowed  to  vanish,  considering  only  positive 
values  of  p0,  the  three  singular  roll  rates  are  identified  as 


Po  =  W1 

p0  =  (3wj  +  w2)/4  (4.12) 

Po  =  («i  +  w2>/2 

which  are  identical  to  the  singular  roll  rates  resulting  from  the  case  of  constrained  rolling  motion  of  Chapter  III. 
Since  the  functional  forms  of  the  induced  yawing  moment,  and  induced  roll  moment  are  quite  similar,  this  result 
is  not  surprising. 


B.  RESONANT  SOLUTIONS 


When  the  rolling  velocity  is  in  the  vicinity  of  the  above  singular  roll  rates,  additional  resonant  solutions  must 
be  generated  which  are  valid  in  the  vicinity  of  these  points.  If  Equation  (4.6)  is  considered  in  the  vicinity  of  funda¬ 
mental  resonance,  and  if  the  0j  are  interpreted  as  cjjT0,  it  is  evident  that  the  term  0]K4  exp[4i(0,  -0)]  will  be 
stationary  for  perfect  resonance,  p  =  W],  For  roll  rates  in  the  vicinity  of  resonance,  p0  ^  to,,  this  term  will,  at 
most,  be  a  slowly  varying  function  of  time,  T2,  and  hence  will  give  rise  to  secular  terms  in  the  expansion  p(t). 
Therefore,  in  order  that  the  solution  remain  valid  for  large  time,  it  is  required  that 


dPo 

dT2 


a2p0  +  °o  +  °lkl  Sin  4(0 1  -0)  =  0 


(4.13) 


9 


Equation  (4.6)  then  becomes 


3p2  9p2  3P2  _  f  4  i4(*2-<M  3  i<3*.+; 

2  Po  ~  a'lm[  2  K‘  26 


-40) 


+  6K2K2e2i,J'*-^-2^  +  (4.14) 


which  displays  a  particular  solution  valid  at  p0  =  to. 


°lk2 


4o,kjk2 


p2(t)  =  4^0 -«2)  COS4(^2  '♦>  +  4p0  -"3w"j"-  cj2  C°S(3*>  +  4/2  ~40) 


3glklk2 
2p0  “  to |  -  co2 


Cos  2(0  ]  +  02  -  20)  + 


40,k,k2 


4p0  to,  3to2 


2^-  Cos  (0, +302 -40)  (4.1  S) 


In  like  manner,  approximate  solutions  which  are  valid  at  zero  spin,  p0  =  to,  +  co2,  and  one-half  resonant  spin, 
p0  =  (3co,  +  co2)/ 4  may  also  be  generated.  For  zero  spin,  the  variation  of  the  first-order  solution  p0(t)  is  given  by 


-jp  =  <*2Po  +  °0  +  aiklk2  Sin  2($i  +  02  _  2 0) 


(4.16) 


while  the  second-order  term  in  the  expansion  is  found  to  be 


°lk'l 


°lk2 


p2(t)  =  77 - 7  Cos  4(0.  -0)  +  - r  Cos  4(02 

2  4(p0-.w1)  1  '  4(p0-to2)  KV2 


-0) 


4°t  kik2  „  «.  4aik]k2 

%p0-3cu,-to2  Cos  (30,  +02-40)  +  ^  ^  ~ ~  Cos  (0,  +  302  -  40)  (4.17) 


Finally,  the  first-order  rolling  motion  at  one-half  resonant  spin  is  governed  by 

^Po  2  a  ^ 

■jp  =  °2Po  +  a0  +  4<7,k,k2  Sin  (30,  +  02  -40) 


(4.18) 


and  the  second-order  solution  becomes 


alki 


fflk2 


p2(t)  =  4fp~0  -to,)  C°S4(^  +  4(p0  ~)  C0S4W2  "♦) 


3k(k 


2c? 


1*2 


_  A  A 


4k,  k^ 


+  2p0-co,  -«2  Cos  2(0,  +  02  ~  20)  +  4po  _  —  -3W2  Cos(0,+302  -40)  (4.19) 

It  should  be  emphasized  that  Equations  (4.13)  through  (4.19)  are  valid  for  the  discrete  roll  rates  only. 


24 


A  direct  comparison  of  the  approximate  solutions  to  Equation  (4.1)  using  the  method  of  multiple  scales  with 
direct  numerical  solutions  is  complicated  by  the  singular  nature  of  the  problem  at  the  three  discrete  roll  rates  men¬ 
tioned  above.  Typical  roll  rates  of  interest  for  operational  air-launched  ordnance  range  from  zero  spin  to  the  onset 
of  Magnus  instability  which  is  usually  several  times  resonance.  Therefore  a  comparison  of  the  approximate  solutions 
of  this  chapter  with  direct  numerical  integrations  of  Equation  (4.1),  over  spin  rates  of  interest  will  necessitate  the 
use  of  multiple  resonant  and  nonresonant  solutions  which  must  be  “matched”  in  an  appropriate  manner  such  that 
continuous  coverage  of  the  roll  spectrum  is  achieved. 


In  practice,  this  matching  is  achieved  through  the  introduction  of  a  detuning  parameter  X,  which  is  used  to 
slightly  detune  the  singular  solution  as  the  roll  rate  varies  about  the  singular  roll  rate.  For  example,  in  the  case  of 
fundamental  resonance,  the  roll  rate  is  detuned  according  to 

p  =  coj  +  Xe2  (4.20) 


where  v  =  0(1).  Substitution  of  Equation  (4.20)  into  Equation  (4.13)  yields 


dpO 

dT, 


°2P0  +  °0  +  °!  kl  Sin  —  00  —  XT2) 


(4.21) 


Thus,  through  the  introduction  of  the  parameter  X,  the  “matched"  asymptotic  expansions  of  Equation  (4.5)  can  be 
used  to  describe  the  nonlinear  rolling  motion  over  the  entire  spectrum  of  roll  frequencies. 

Figures  12  through  18  display  the  missile  induced  rolling  motion  for  yawing  motion  of  various  amplitudes. 
Missile  physical  properties  and  aerodynamic  characteristics  used  in  the  calculations  were  shown  in  Table  1 .  For  the 
purposes  of  these  calculations,  approximate  solutions  generated  by  the  method  of  multiple  scales,  were  matched 
at  points  which  were  chosen  to  lie  midway  between  resonant  roll  rates.  While  more  elaborate  methods  may  be 
devised,  the  present  method  proved  reasonably  successful  as  evidenced  by  the  agreement  displayed  in  Figures  12 
through  18  for  a  variety  of  initial  conditions  and  design  roll  rates.  These  results  display  the  dependence  of  the 
roll  rate  history  on  initial  conditions  and  angle  of  attack  for  a  passage  through  resonance  to  the  design  roll  rate 
of  3cj[.  As  the  amplitude  of  the  yawing  motion  increases,  the  rolling  motion  becomes  dominated  by  the  large 
induced  roll  effects,  and  the  validity  of  the  perturbation  analysis  becomes  questionable  at  angles  of  attack 
in  the  neighborhood  of  30  degrees.  In  addition,  the  inaccuracy  of  the  asymptotic  expansion  of  Equation (4.5)  for 
large  time,  is  evidenced  by  the  growing  phase  difference  oetween  the  approximate  and  exact  solutions. 

Missile  roll  behavior  ip  the  vicinity  of  resonance  may  be  demonstrated  by  an  inspection  of  Equation  (4.13) 
together  with  an  analysis  of  the  following  numerical  results.  Figures  19  through  21  show  missile  roll  rate  histories 
passing  from  resonance  (p0  =  wj )  to  a  design  roll  rate  of  3u>j ,  under  the  influence  of  progressively  larger  yawing 
motion  which  is  predominantly  precessional.  As  predicted  by  Equation  (65),  the  roll  rate  approaches  its  design 
value  with  only  minor  variations  even  as  the  amplitude  of  the  yawing  motion  approaches  some  30  degrees. 


Stability  of  the  rolling  motion  at  resonance  is  demonstrated  by  the  results  of  Figures  22  through  27  where 
the  missile  design  roll  rate  was  chosen  as  pss  =  .  The  diminished  influence  of  large  amplitude  precessional  motion 

is  again  displayed  by  the  results  of  Figures  22  and  23  while  Figures  25,  26,  and  27  indicate  the  drastic  change  in 
the  character  of  the  rolling  motion  under  the  influence  of  large-amplitude  nutational  motion. 
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Figure  16.  -  Missile  Induced  Rolling  Motion  -  Passage  Through  Resonance 
(k,  =0.21,  k2  =0.02,  ar=  13.2°) 
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Figure  1 7.  Missile  Induced  Rolling  Motion  -  Passage  Through  Resonance 
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Figure  18.  -  Missile  Induced  Rolling  Motion  -  Passage  Through  Resonance 
(k,  =  0.41.  k2=  0.05,  oT  =  26.4°) 
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Figure  19.  -  Missile  Roll  Rate  History  (k,  =0.1.  k-%  =  0.2,  aT  =  13.2°) 
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Figure  20.  -  Missile  Roll  Rate  History  (k,  =  0.1 .  k2  =  0.3,  aT  =  22.9°) 
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Figure  21.  -  Missile  Roll  Rate  History  (k,  =  0.10.  k,,  =  0.4.  =  28.6“) 
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Figure  26.  -  Missile  Roll  Rate  History  (k,  =  0.21 ,  k2  =  0.02,  aT  =  1 3.2°) 
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Figure  27.  -  Missile  Roll  Rate  History  (k,  =  0.31.  k,  =  0.03.  a-r  =  19.5°) 


CHAPTER  V 


COMBINED  ROLL-YAW  INTERACTION 


The  results  of  the  preceding  chapters  have  been  instrumental  in  providing  quantitative  insight  into  the 
mechanisms  of  resonance  instability  as  related  to  the  effects  of  roll  orientation-dependent  aerodynamics.  In  addi¬ 
tion,  stability  criteria  useful  to  the  designer  have  surfaced  as  a  by-product  of  approximate  solutions  generated  by 
the  perturbation  methods.  While  the  two  specific  conditions  studied,  i.e.,  free  yaw  with  prescribed  rolling 
motion  and  induced  rolling  motion  with  prescribed  yawing  motion,  have  been  instructive,  neither  can  be  considered 
totally  representative  of  free-flight  conditions.  In  this  chapter,  approximate  solutions  are  sought  which  are  represen¬ 
tative  of  the  full  nonlinear  problem  of  combined  roll-yaw  interaction. 

The  simultaneous  differential  equations  governing  the  combined  roll-yaw  problem  may  be  taken  directly  from 
Chapter  III  (Equation  (3.1))  and  Chapter  IV  (Equation  (4.1)) 

?  -  iP0|  -  <41  =  e2Mi?  +  ie2/i2£  +Xi52e_l47|  (5.1 ) 

and 

p  =  e4(o0  +  o2p)  +  Im(]jV'4*)  (5.2) 

Note  that  the  smallness  parameter  e  has  been  introduced  into  Equation  (5.2)  in  a  slightly  different  manner,  such 
that  terms  on  the  right-hand  side  are  of  equal  order  of  magnitude. 

Approximate  solutions  to  the  system  of  nonlinear  Equations  (5.1)  and  (5.2)  are  again  sought  using  a  generalized 
version  of  the  method  of  multiple  scales  similar  to  that  of  Chapter  IV.  It  is  assumed  that  £  and  p  possess  uniformly 
valid  expansions  of  the  form 


|(t;e)  =  e|i(^l,i^2,T2)  +  e3£3(\|/,,  T2)  +  ... 


p(Ce)  =  p0(T4)+  e4p4(^],  ii/2,0,T4)  +  ... 


(5.3) 


where 

d^ 

—  =Wj,  (5.4) 

and 

Tj  =  e*T0  (5.5) 

The  slow  time  scale  T2  is  chosen  so  as  to  characterize  the  slowly  varying  nature  of  the  modal  amplitudes  while  the 
variable  T4  is  introduced  into  the  roll  expansion  in  order  to  accommodate  the  comparatively  slow  variation  of  spin 
due  to  C6()  In  addition,  the  assumed  expansions  of  Equation  (5.3)  differ  from  previous  expansions  in  that  £  may 
now  exhibit  roll  dependence  through  the  independent  variable  <t>. 
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Equations  (5.3)  may  now  be  substituted  into  Equations  (5.1)  and  (5.2)  to  obtain  partial  differential  equations 
in  the  independent  variables  0 j,  Tj  and  0.  Equation  (5.3)  may  be  incorporated  into  Equation  (5.2),  equating  all 
terms  of  order  e4  to  obtain 


U) 


9p4  dp4  3p4 

1  30t  +  302  +  P°  30 


dPo 

~dT^  +  a°  +  °2Po 


+  0)  lm(|4e"l44>) 


(5.6) 


Equation  (5.6)  is  identical  to  the  roll  partial  differential  equation  which  resulted  from  the  analysis  of  Chapter  IV 
with  the  exception  of  the  fact  that  is  no  longer  constrained  to  be  epicyclic  with  constant  amplitudes  and  fre¬ 
quencies  but  may  vary  in  accordance  with  Equation  (5.1 ). 

Equation  (5.3)  may  be  substituted  into  Equation  (5.1)  and  terms  of  order  e  and  e3  equated  to  obtain  partial 
differential  equations  for  the  solution  of  |[  and  |3 . 

Terms  0(e): 


>*«!)- iPoVdO-Mfo  =  0 


Terms  0(e3): 


^2(?3)-iPo^3)-“0?3  -  -2cu, 


32?1  _  92$i 

30,3T2  2C°2  3023T2 


dtj]  3£j  dw2 

dT2  30j  dT2  90 2 


+iPn 


dT, 


+  Mi^dj)  +  4i2Si  +  Xili3el4',> 


(5.7) 


(5.8) 


where 


'J  =  Ui 


3 

90! 


w2 


902 


(5.9) 


and  the  epicyclic  frequencies  cjj  are  taken  to  be  functions  of  slow  time  T2,  to  characterize  the  slow  variation  of  the 
cjj  with  mach  number,  dynamic  pressure,  etc.  Thus  Equations  (5.6)  through  (5.9)  represent  the  governing  equations 
for  fp  p0,|3,and  p4. 


A.  NONRESONANT  SOLUTIONS 


The  general  solution  to  the  homogeneous  Equation  (5.7)  has  the  usual  epicyclic  form  and 
directly  as 


may  be  written 


where 


and 


2 


l  =  E 

(5.10) 

j=l 

Ki  =  k/fli 

(5.11) 

0j  =  WjT0 

(5.12) 
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Equation  (S.10)  may  be  substituted  into  Equation  (S.8)  to  obtain  the  partial  differential  equation  governing  |j 
resulting  in 


t)KiK1e'(‘,'!*‘‘*11  *  3K1cieK4*'*'':*11]  (5.13) 

where 

dK,  /  \ 

w,  =  i(P0  -  2Wl)  -gj-  +  i  f^co,  +  p2  -  -^ij  K,  (S.14) 

and 

dK,  /  dco2  \ 

*2  =  +  >^1^2  +#12  ~^)K2  (515) 

In  order  to  obtain  a  uniformly  valid  solution  to  Equation  (5. 1 3)  that  remains  bounded  as  T0  -*•  it  is  required  that 
the  secular  terms  vanish  (wj  =W2  =  0),  which  results  in  the  differential  equations  characterizing  the  slowly 
varying  Kj’s 

dK,  /  dco,  \ 

i(P0  -  2w, )  +  i  =°  ($16) 

dK2  /  d(x>2  \ 

i(Po-2w2)-^  +  2  _dT7jK2  =  0  (517) 

With  these  secular  terms  removed,  the  particular  solution  to  Equation  (5.13)  becomes 


where 


(5.18) 


_ -XjK? _ 

(4p0  -  3wj)2  -  P0(4p0  ~  3wj)  +  Wq 

_ -xM _ 

(4p0  -  3w2)2  “  P0(4p0  -  3w2)  + 

_ -3XiKfK2 _ 

(4p0  -  2ojj  -  uj2)  -  P0(4p0  -  2u>j  -  oo2)  +  w l 


(5.19) 


and 


Oa  = 


-3xiK,K| 


(4po  co,  2gj2)  Pq(4Pq  Wi  2w2)  +  cog 
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Returning  to  the  differential  equation  of  rolling  motion.  Equation  (5.10)  may  be  substituted  into  Equation  (5.6), 
and  again  eliminating  the  secular  term,  the  slow  variation  of  p0(T4)  is  given  by 

dp0 

~dT~  +  °0  +  °2P0  =  0  (5-20) 

U  14 

The  particular  solution  to  Equation  (5.6),  with  the  secular  term  removed,  may  then  be  written 
P4(0  =  Cos  4  (i£j  -  0)  +  02  Cos  4  (02  ~  <t>)  +  03  Cos  (3^j  +  1 p2  -  4 <p) 

+04  Cos  2  (0j  +  -  20)  +  05  Cos  (0j  +  302  -40)  (5.21) 


with 


-  gtki 

^  4(p0-w,) 

°lk2 

4(Po  -  0)2) 
4o.kfk2 

Q  —  _ 1  1  * 

^  *Po  ~  3cOj  ~  co 2 
3o,k^k2 

/3  —  _ 1  1  * 

4  2p0  -  cjj  -  w2 

4o,  k.k2 

0,  =  - - i — 1— : - 

4p0  -  u>]  -  3w, 

and 

0j  =  WjT0  +  flj 


(5.22) 


(5.23) 


Inspection  of  Equations  (5.10)  through  (5.23),  which  constitute  the  first  and  second  approximate  solutionsfor  p(t) 
and  j-(t),  reveals  that  the  full  problem  of  combined  roll-yaw  interaction  may  be  treated  as  a  superposition  of  the 
two  previously  treated  problems,  so  long  as  the  assumptions  regarding  the  perturbation  nature  of  the  problem 
remain  valid. 


B.  RESONANT  SOLUTIONS 

Similar  to  the  behavior  investigated  in  Chapters  III  and  IV,  the  higher-order  approximate  solutions  for  coupled 
yawing  and  rolling  motion,  (|3,  p4),  become  singular  as  the  roll  rate  approaches  the  distinct  values  of  zero  spin, 
resonance  (p  =  Wj)  and  approximately  one-half  resonant  spin.  This  difficulty  is  evidenced  by  the  singular  behavior 
of  the  terms  of  Equations  (5. 19)  and  (5 .22)  as  p0  approaches  one  of  these  critical  roll  rates.  Further  evidence  may  be 
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obtained  through  an  observation  of  the  nonhomogeneous  terms  of  Equations  (5.6)  and  (5.13)  as  the  roll  rate 
approaches  the  above  critical  values.  Selected  high  frequency  terms  become  stationary  or,  at  best,  are  slowly  varying 
functions  of  time  (T2,T4),  which  give  rise  to  secular  terms  in  the  expansions  for  p  and  |.  Perturbation  solutions 
which  are  valid  in  the  vicinity  of  these  singular  roll  rates  must  be  generated  individually  and  these  solutions 
“matched"  to  obtain  continuous  approximate  representation  of  the  motion.  Solutions  valid  in  the  vicinity  of 
resonance  are  generated  in  what  follows.  For  the  sake  of  brevity,  derivation  of  solutions  valid  at  p0  =  0, 
(3cjj  +  a>2)/4  are  omitted  but  may  be  easily  derived  in  a  manner  similar  to  the  case  of  resonant  spin,  p  =  go,  . 

At  resonance,  the  0j  and  0  are  interpreted  as  cojTq  and  PqTq  respectively  and  Equation  (5.13)  may  be 
rewritten,  grouping  stationary  terms,  as 


V  2(?3)  -  iP0y(?3)  -  cjq|3  =  [wj  +  x,K3ei40<>]ei0'  +w2e^ 

+  Xil^e*40'3^  +  3KjiC +  skfoe*40-20'-^]  (5.24) 
In  order  for  the  solution  to  Equation  (5  24)  to  remain  uniformly  valid,  it  is  required  that  the  secular  terms  vanish 

w,  +  xi  Kje1400  =  0  (5.25) 

w2  =  0  (5.26) 

The  solution  of  Equation  (5.24),  with  secular  terms  removed,  may  now  be  obtained  as: 

?3(t)  =  a2ei(40-30>>  +  aje1'40-20*"0^  +  a4ei,40'0>"2^)  (5.27) 

Thus  the  resulting  asymptotic  expansion  for  ?(t),  valid  at  fundamental  resonance,  is  written  as 

|(t;e)  =  ei1(^1,^2,T2)  +  e3I3(^1,^2,0.T4)  (5.28) 

with  and  f3  given  by  Equations  (5.10)  and  (5.27)  respectively. 

The  solution  to  the  nonlinear  rolling  motion  at  resonance  is  generated  by  a  consideration  of  Equation  (5.6) 
as  p  -*  cj[ .  Grouping  stationary  (slowly  varying)  terms,  Equation  (5.6)  may  be  expanded  and  rewritten  as: 


dp4  9p4  9p4  dPo 

W1  ^  +  +  Po  dl7  +  Po  +  °2P0  +  °lkl  Sin  4  (0 ,  -  0O ) 


+a^k2  Sin  4  (02  -  0)  +  4k,  k2  Sin  (30,  +  02  ~  40) 

+6k3k2  Sin  2  (0]  +  02  -  20)  +  4k,  k2  Sin  (0,  +  302  -  40)J  (5.2C>) 


Elimination  of  the  slowly-varying  secular  terms  yields  the  governing  equation  for  p0: 

dPo  , 

=  a0+  a2p0  +  o,k|  Sin  4(0,  -0O) 


v-  . 


>{**'*■'  ‘X-< 


The  uniformly  valid  solution  for  p4  is  then  solved  for  directly  and  is  found  to  be: 


P4(0  =  02  Cos  4  ($2  ~  0)  +  03  Cos  (30 1  +  02  “  40> 

+04  Cos  2  (0,  +  02  -  20)  +  0S  Cos  (0,  +  302  ~  40) 


(5.31) 


C.  STABILITY  AT  RESONANCE 

Equations  (5.27)  through  (5.31)  comprise  approximate  solutions  to  the  combined  roll-yaw  motion  valid  in 
the  vicinity  of  resonance.  While  these  results  contain  considerable  information  regarding  the  character  of  the 
coupled  motion  at  resonance,  they  are  nonlinear,  time-varying,  coupled  differential  equations  that  yield  little  quanti¬ 
tative  information  without  recourse  to  direct  numerical  integration.  Additional  insight  into  the  stability  of  the 
motion,  in  close  proximity  to  resonant  spin  may  be  gained  through  a  small  perturbation  analysis  of  Equations  (5.27) 
through  (5.31). 

Equation  (5.25)  may  be  expanded  into  real  and  imaginary  components  to  yield  the  differential  equations 
governing  the  slowly  varying  nutation  amplitude  and  phase 


dk, 

dfj 


/Mi^l  +  H2  ~ 

\  P0  2wl 


-)k» 


Xi^i 

P0  “  2<o, 


Sin  4(0-0]) 


(5.32) 


dp0  Xjki 

Wi'TTM0*4*'*1'  (533) 

The  governing  differential  equation  of  rolling  motion  at  resonance  is  taken  from  Equation  (5.30)  as 

=  a0  +  a2p0  +  a,k?  Sin 4(0-0,)  (5.34) 

An  investigation  of  the  combined  rolling  and  yawing  motion  in  the  neighborhood  of  p0  =  w,  proceeds  by  consider¬ 
ing  the  phase  0-0,  to  be  essentially  stationary,  or  at  most  a  slowly  varying  function  of  T2  at  resonance.  Allowable 
perturbations  about  resonance  are  made  possible  through  the  introduction  of  the  detuning  parameter  X,  such  that 


p0  =  w,  +  e2\ 


(5.35) 


where  X  is  of  order  unity.  Substitution  of  Equation  (5.35)  into  Equations  (5.32)  through  (5.34),  with  the  interpre¬ 
tation  of  0,  as  w,T0, yields  the  following  nonlinear  system 


dkl  _  /M i«i  +  H2  ~  «?\  _  Xikj 

dT2  ~  \  p0  -  2W,  )  "  p0  -  2w, 


(5.36) 


dr  _  „  4Xi 
dT2  "  P0  -  2co, 


k?  cot  r 


(5.37) 


^  -  °lkl  Sinr 


(5.38) 


40 


where 


Figure  31.-  Root  Locus  at  Resonance,  k. 


only  mild  instabilities  associated  with  these  roots  with  increasing  frequency  as  the  nutation  amplitude  increases. 
However,  strongly  divergent  instabilities  are  possible  due  to  a  third  real  root,  which  may  be  large  and  positive. 

The  analytical  results  of  Equation  (5A6)  were  verified  by  direct  numerical  solutions  of  Equations  (5.1)  and 
(5.2)  with  the  results  shown  graphically  in  Figures  33  through  38.  Mild  instabilities  about  resonance  are  depicted 
by  the  total  angle  of  attack  and  roll  rate  histories  in  Figures  33  through  35  for  various  values  of  nutation  amplitude. 
As  the  amplitude  of  the  motion  is  increased,  the  mild  oscillatory  behavior  of  the  nutation  arm  becomes  more 
pronounced  in  agreement  with  the  perturbation  analysis. 

The  oscillatory  behavior  of  the  nutation  amplitude  associated  with  the  solutions  of  Figures  33  through  35  is 
shown  in  Figure  36.  Frequencies  characteristic  of  this  behavior  are  seen  to  be  in  agreement  with  the  root  loci  of 
Figures  28  through  30.  Note  that  the  missile  roll  rate  has  departed  significantly  from  its  equilibrium  value,  p  =  coj, 
where  the  perturbation  assumption  is  no  longer  valid.  In  Figure  37,  the  missile  stability  in  the  vicinity  of  resonant 
spin  is  shown  to  be  quite  good  for  initial  conditions  producing  small  amplitude  motion.  In  contrast,  Figure  38 
illustrates  the  sensitivity  of  the  missile  stability  to  small  changes  in  initial  conditions  for  large  amplitudes  motion. 
Note  that  although  both  solutions  of  Figure  38  have  identical  initial  conditions,  with  the  exception  of  small 
variations  in  ,  their  angular  motion  and  resultant  spin  rates  are  significantly  different. 

While  not  directly  applicable  to  the  small  perturbation  results,  Figures  39  and  40  show  total  angle  of  attack 
and  roll  rate  histories  for  configurations  with  fin  cant  with  a  resulting  design  roll  rate  of  3coj.  The  results  of 
Figure  39  show  stable  missile  angular  motion  when  initial  conditions  are  chosen  in  accordance  with  previous  root 
loci  predicting  stable  motion  (kj  =  0.3,  Tq  =  3.6).  The  steady-state  roll  rate  is  seen  to  exponentially  approach  its 
design  value  with  apparently  little  influence  from  the  missile  angular  motion.  However,  when  the  initial  phase  angle 
r0  is  chosen  such  that  the  dominant  real  root  of  the  linear  perturbed  system  is  large  and  positive,  the  divergent 
angular  motion  of  Figure  40  is  obtained  together  with  a  rapid  divergence  of  the  roll  rate  to  a  value  well  beyond 
the  intended  design  value  of  pss  =  3cj1  .  The  critical  value  of  spin  producing  the  onset  of  magnus  instability  for  this 
particular  configuration  was  calculated  to  be  (p/coi)cm  =  11.0.  The  maximum  roll  rate  produced  by  the  resonance 
instability  was  approximately (p/ ( )|max  9.2,  justbelow  that  spin  needed  for  magnus  instability. 
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Figure  34.  -  Missile  Stability  at  Resonance 
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Figure  35.  -  Missile  Stability  at  Resonance 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 


The  free-flight  angular  motion  of  a  symmetric  cruciform  missile,  in  the  presence  of  nonlinear-roll-orientation- 
dependent  aerodynamics  has  been  treated  in  a  manner  amenable  to  deriving  practical  stability  criteria.  An  aero¬ 
dynamic  force  and  moment  system,  consistent  with  arguments  of  rotational  and  reflectional  symmetry,  have  been 
incorporated  into  the  equations  of  motion.  The  method  of  multiple  scales  has  been  employed  to  obtain  approxi¬ 
mate  expressions  prescribing  the  missiles  angular  motion  for  the  following  conditions: 

1 .  Free-flight  yawing  motion  in  the  presence  of  nonlinear-roll-orientation-dependent  aerodynamics  with 
prescribed  rolling  motion 

2.  Free-flight  rolling  motion  in  the  presence  of  linear  damping  and  nonlinear  induced  roll  moments 

3.  Aerodynamically  coupled,  free-flight  yawing  and  rolling  motion  both  away  from,  and  in  the  neighborhood 
of  resonant  spin. 

For  the  case  of  free-flight  yawing  motion  with  a  priori  rolling  motion,  an  application  of  the  method  of  multiple 
scales  yields  first-approximate  solutions  identical  to  those  generated  by  the  method  of  averaging.  Critical  roll  rates 
are  identified,  and  solutions  valid  in  the  vicinity  of  resonance  are  derived.  Approximate  solutions  show  good 
agreement  with  direct  numerical  solutions  for  moderate  angles  of  attach  and  as  such  should  provide  useful  stability 
criteria. 

Nonlinear  rolling  motion  with  prescribed  yawing  motion  was  analyzed  using  a  generalized  version  of  the 
method  of  multiple  scales.  Singular  roll  rates  of  zero-spin,  resonance,  and  one-half  resonant  spin  were  identified 
and  solutions  valid  in  the  neighborhood  of  these  roll  rates  derived.  Comparisons  with  direct  numerical  integrations 
were  accomplished  by  matching  resonant  and  nonresonant  solutions  throughout  the  roll  rate  spectrum.  Approxi¬ 
mate  solutions  were  shown  to  remain  valid  for  moderate  angles  of  attack  below  30  degrees. 

The  complete  problem  of  coupled  nonlinear  rolling  and  yawing  motion  was  investigated  and  approximate 
solutions  describing  the  motion  were  found  to  be  simultaneous  solutions  of  the  prescribed  roll  and  yaw  cases. 
Solutions  valid  in  the  vicinity  of  resonance  were  derived  and  combined  with  small  perturbation  methods  to  derive 
useful  stability  criteria. 

The  above  analyses  have  clearly  shown  the  existence  of  resonance  instabilities  due  solely  to  the  inherent  lack 
of  complete  axial  symmetry  introduced  by  the  fin  cruciform.  Through  the  use  of  the  derived  approximate  solutions, 
stability  criteria  may  be  used  to  investigate  the  frequency  of  occurrence  or  susceptibility  of  candidate  configurations 
to  resonance  instabilities. 
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APPENDIX 


MAPLE-SYNGE  THEORY 


The  essence  of  the  Maple -Synge  theory  lies  in  its  representation  of  the  fluid  forces  and  moments  as  power 
series  expansions  in  the  cross  velocity  components  and  in  its  arguments  of  the  symmetry  of  these  forces  and 
moments  under  coordinate  transformation. 

For  purposes  rf  illustration,  the  derivation  of  the  static  forces  and  moments  are  considered  here.  For  a  more 
complete  treatment,  the  reader  is  referred  to  Reference  12. 

The  complex  transverse  force  and  moment  coefficient  expansions,  by  the  “aerodynamic  hypothesis,"  may  be 
written  in  the  missile-flxed  system  in  the  complex  angle  of  attack  1-,  where 


as 


v  +  iw 
V 


(A-l) 


and 


(Cy  +  iCz)  =  £  Cyt1? 
hi 


(A-2) 


(Cm+iCn)  =  £  <*«$*!*  (A-3) 

M 

The  scalar  V  is  the  magnitude  of  the  total  velocity  and  (“)  denotes 
complex  conjugation.  The  coefficients  Cjj,  C^e,  are,  in  general,  complex 
and  may  also  be  functions  of  spin,  Mach  number  and  Reynolds  number. 

The  consequences  of  tetragonal  symmetry  may  be  realized  by  a 
consideration  of  the  cross  force  and  moment  coefficient  expansions  in  a 
new  coordinate  system  (x,  yj ,  zf)  obtained  by  a  90°  rotation  about  the 
x-axis.  The  rotational  symmetry  of  the  missile  requires  that  the  force 
and  moment  coefficients,  as  viewed  in  both  systems  must  be  related  by 

(Cy]  +iCZ[)  =i(Cy+iCz) 

+'C„i)  =  i(Cm  +  iCn)  (A-5) 

with 

Si  =  iS  (A-6) 

and 

I,  =  -il  (A-7) 

A-l 


Now  assuming  that  the  coefficients  Cjj  and  Ckt  remain  invariant,  Equations  (A-2)  and  (A-3)  may  be  substituted 
into  Equations  (A4)  and  (A-5)  to  obtain 


22  CyU^ttjV  =  i  22  (A-8) 

i.j  i,j 

22cMtt1)k(t1)*-i£ck8|kf!  (A-9) 

M  k,s 

Since  the  power  series  of  Equations  (A-8)  and  (A-9)  must  be  equal  term-by-term,  consideration  of  Equations  (A-6) 
and  (A-7)  requires  that 


iH-l  =  i 


and 


:k-c-l 


(A-IO) 


(A-ll) 


or  equivalently 


i-j- 1  =  k-8-1  =  4N 


(A-12) 


where  N  is  an  integer. 


Therefore,  the  tetragonal  symmetry  requires  that  only  those  values  of  i,  j,  k,  and  8,  in  compliance  with  Equa¬ 
tion  (A-12),  may  occur  in  the  cross  force  and  moment  expansions. 


*2 


The  consequences  of  reflectional  symmetry  may  be  deduced  by  a 
consideration  of  the  form  of  Equations  (A-2)  and  (A-3)  in  a  coordinate 
system,  (x,  y2,  z2)  obtained  by  reflection  about  the  x-y  plane.  However, 
this  reflection  has  changed  the  right-handed  coordinate  system  (x,  y,  z)  to 
a  left-handed  (x,  y2,  z2)  system.  Thus,  the  moment  coefficient  (Cm2  +  iC„2) 
must  be  redefined  in  the  (x,  y2,  z2)  system.  Maple  and  Synge  (A-2)  adopt 
the  definition  that  a  rotation  is  positive  if  it  corresponds  to  a  cyclical 
rotation  of  the  axes. 

With  this  in  mind,  the  force  and  moment  expressions  in  the  two 
systems  are  related  by 


(Cy2  +  iC^)  =  (Cy  +  iCz)  (A- 13) 


and 


(Cm2+iC„2)  =  -(Cm+iCn) 


(A- 14) 


with 


$2  =  $ 


(A-15) 


A-2 


1 


Substitution  of  Equations  (A-2)  and  (A-3)  into  Equations  (A-l  3)  and  (A-14)  yields 


£  Cytt2)V«2 )j  =  £  Cjjft)^)1  (A-l 6) 

i.j  i.j 


L  c™kMCh)v  =  -£  wdV 

k,8  k,C 


Incorporation  of  Equation  (A-l  5)  leads  to  the  conclusions 


and 


C  =  -C 

mk$  mkC 


(A-l  7) 


(A-l  8) 


(A-l  9) 


Thus,  reflectional  symmetry  requires  the  force  coefficients  Cjj  to  be  real  and  the  moment  coefficients,  CmkC  to  be 
pure  imaginaries. 

In  its  most  general  form,  then,  the  force  and  moment  coefficient  expansions  for  a  symmetric  cruciform 
missile  may  be  written 


(Cy  +  iCz)  =  £  CNj(£)j  +  1+4N(F)j  (A-20) 

N,  j 

(Cm  +ic„)  =  2]  CmNe(t)e+1+4N(«B  (A-2 1 ) 

N,8 

Expansions  of  Equations  (A-20)  and  (A-2 1 )  for  various  values  of  N  exhibit  the  form  of'the  roll-orientation- 
dependent  aerodynamics.  For  the  case  N  =0,  Equations  (A-20)  and  ( A-2 1 )  yield 


(Cy  +  iCz)^  =  2]  C0.62j$  (A-22) 


(Cm+iCn)o  =  £  Cm0iS6Jet  (A-23) 

c 

where  6  2  =  |f.  Since  all  the  terms  in  the  expansions  of  Equations  (A-22)  and  (A-23)  are  collinear  with  £,  it  is  clear 
that  the  N  -  0  terms  are  incapable  of  discerning  the  presence  of  the  fins.  However,  consideration  of  the  N  =  ±1 
terms  may  be  used  to  generate  the  roll  orientation-dependent  terms  harmonic  in  Ay.  For  N  =  ±1  the  force  and 
moment  expansion  becomes* 


'Here  the  force  coefficients  (C’s)  have  been  replaced  by  the  more  conventional  notation  (C^’s)  with  the  subscript  N  denoting 
normal  force. 


LCn...'*Cn,.|!,2*CnI.284'‘'' 

¥ 

(A-24) 

<c«  ♦iCn)N,11  ■  s«4  ♦  • 

-]PH 

[C",^C»,,62*C".,,84*- 

(A-25) 

Equations  (A-24)  and  (A-25)  may  be  put  in  standard  form  if  the  force  and  moment  vectors  are  transformed 
to  the  nonrolling  system  through  the  transformations'. 

(C^  +  iC?)  =  e^Cy  +  iCz)  (A-26) 

(C~+iCli)  =  e^(Cm+iCn)  (A-27) 

and 

I  =  ei0$  (A-28) 

with  the  result 

(C?  +  iC*)N=±I  =  [CN_1i3  4  Cn_1>3  4  Cn_i>452  +  CN  ]_s54  +  . .  js2e-i47| 

+[CN(  o  +  CNi>i52  +  CN)/  +  . .  .]64ei47F  (A-29) 

^  +iCn)N=tl  =  K.,3  +C-n-,452  +C-n-,s64  +  '^^1 

+  [Cmt,o  +  C»,..62  +  cn.,.2«4  +  -]«4ei41'?  (A-30) 

where 

y  =  e-<p  (A-31) 

Inspection  of  Equations  (A-29)  and  (A-30)  clearly  reveals  that  the  lowest-order  terms  capable  of  exhibiting  roll 
orientation-dependence  are  third  order  in  5. 

Combining  the  expansions  for  N  =  0,  1,-1,  static  aerodynamic  force,  and  moment  coefficients  to  third  order 
in  8  may  be  obtained 

(C~y  +  iCj)  =  [cNq  q  +  C^S2  +  CN_]  382e"i47  +  . .  ]}  (A-32) 

and 

(Cffi  +iC~n)  =  [crao  o  +Cmo  t82  352e-i4"  +  . .  .]*  (A-33) 


Since  all  the  Cmjj’s  are  pure  imaginary,  Equations  (A-32)  and  (A-33)  may  be  written  in  the  more  conventional 
aerodynamic  notation: 

(CT+K~)=[cN  +Cnj82+Cn  282e-i4l'+...]|  (A-34) 

L  a  ” 

(Cm  4 iC75)  =  +  Cm^2 5 2  +  C^V*4’  +...]?  (A-35) 

A- 4 


Static  aerodynamic  expansions  associated  with  the  missile  rolling  motion  may  be  derived  using  the  same 
arguments  of  symmetry  as  in  the  case  of  missile  transverse  aerodynamics.  The  static  roll  moment  coefficient,  Cc,  is 
assumed  to  have  a  similar  power  series  representation: 

C*  =  E  CU^J  (A’36) 

ij 

where  the  Cy  are  real.  Following  Maple-Synge  symmetry  arguments,  if  Ct  is  viewed  in  the  coordination  system  of 
Figure  A-l ,  Cc  must  remain  invariant,  or 

Cc  =  CCj.  (A-37) 

Equation  (A-37)  may  be  expanded  through  the  use  of  Equations  (A-6),  (A-7),  and  (A-36)  to  obtain 

E  Ckc(S)k(?)C  =  E  Cks(?i)k(?|)e  =  E  Ck6(i)k'c($)k(?)6.  (A-38) 

ke  kt  ke 

Since  the  two  series  must  be  equal  term  by  term,  it  is  required  that 

ik-s  -  1  (A-39) 


or  equivalently 


k  -  8  =  4N  N  =  0,  +1,  ±2,  (A-40) 

Thus  the  roll  moment  coefficient  expansion  must  take  the  form 


C6  =  E  CNc(l)S+4N(?)e  N  =  0,  ±1,  ±2,  .... 

N,  5 

Now  Equation  (A-41)  may  be  rewritten  in  the  slightly  different  form 


where  6 2  =  £f.  Now  the  terms  with  index  N  =  0  may  be  extracted  from  the  summation  to  yield 


c*  =  E  co +  E  E  s2?[cn^4N  +  c.NEr4N] . 

5=0  N  "1  6=0 


(A-41) 


(A-42) 


(A-43) 


which  is  a  statement  of  the  linear  superposition  of  the  static  roll  driving  terms,  which  are  functions  of  yaw  amplitude 
only,  and  the  induced  roll  terms  which  are  functions  of  (•  and  its  conjugate. 

Expressing  the  complex  angle  of  attack  %  in  the  missile  fixed  system  as  Sen,  the  roll  moment  coefficient 
expansion  may  be  written 


Cs  =  E  C06«2t  +  E  54N  E  S2S{Cwe>4Ny  +  C_N6e-i4N^}  (A-44) 

2=0  N  =  1  2  =  0 
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n 


Since  C,  must  remain  a  real  quantity  the  condition 

Cc  =  C,  (A-45) 

may  be  imposed  and  the  two  resulting  infinite  series  equated  term  by  term  with  the  results 

C0,  =  C_0c  (A-46) 

and 

^N2  =  C_NC  (A-47) 

Incorporating  Equations  (A-46)  and  (A-47)  into  Equation  (A-44)  letting  CNf  =  DNcei4s,  the  roU  moment  coefficient 
expansion  yields: 


=  E  co^2t  +  E  * 

N=1 


2=0 


4N 


2E  &2cDn,  Cos  4(/3  +  N7) 
L  2=0 


(A-48) 


where  0  is  a  constant  angle  related  to  the  orientation  of  the  fin  cruciform  in  the  missile-fixed  coordinate  system.  If 
the  cruciform  is  oriented  coincident  with  the  y-z  axes  (as  is  often  the  case)  symmetry  considerations  require  that 
the  induced  roll  moment  vanish  at  7  =  0,  n/2,  etc.,  therefore 

(3  =  n/8  (A-49) 

Thus  the  final  form  of  the  Maple-Synge  static  roll  moment  coefficient  expansion  takes  the  form 


00  00  00 

C2  =  E  C0  5526  +  ^  64N  ^  62EBfjc  Sin  4Ny  (A-50) 

2=0  N=1  2=0 


The  lowest-order  terms  in  the  expansion  may  be  generated  by  a  consideration  of  the  terms  associated  with 
C  =  0, 1  to  obtain 

Cc  =  Cqo  +  C01S2  +  ...  +  B1054  Sin  47  +  Bu56  Sin  4y  +  B2068  Sin  87  +  B2i610  Sin  87  +  ...  (A-51) 

or  written  in  more  conventional  aerodynamic  notation 

Ce  =  Ce  +  C,  52  +  ...  +  Cp  .S4  Sin  4y  +  Cc  ,S6  Sin  4y 
*  0  V  “76  7»6 

+C„  .5®  Sin  87  +  Cc  ..610  Sin  87  +  ...  (A-52) 

Ky  «8  v7«10 
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